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Abstract. The shell-model Monte Carlo (SMMC) technique transforms the traditional 
nuclear shell-model problem into a path-integral over auxiliary fields. We describe 
below the method and its applications to four physics issues: calculations of sd-pf- 
shell nuclei, a discussion of electron-capture rates in p/-shell nuclei, exploration of 
pairing correlations in unstable nuclei, and level densities in rare earth systems. 



I INTRODUCTION 

Studies of nuclei far from stability have long been a goal of nuclear science. Nuclei 
on either side of the stability region, either neutron- rich or deficient, are being 
produced at new radioactive beam facilities across the world. At these facilities, 
and with the help of advances in nuclear many-body theory, the community will 
address many of the key physics issues including: mapping of the neutron and 
proton drip lines, thus exploring the limits of stability; understanding effects of the 
continuum on weakly bound nuclear systems; understanding the nature of shell gap 
modifications in very neutron-rich systems; determining nuclear properties needed 
for astrophysics; investigating deformation, spin, and pairing properties of systems 
far from stability; and analyzing microscopically unusual shapes in unstable nuclei. 

The range and diversity of nuclear behavior, as indicated in the above list of on- 
going and planned experimental investigations, have naturally engendered a host of 
theoretical models. Short of a complete solution to the many-nucleon problem, the 
interacting shell model is widely regarded as the most broadly capable description 
of low-energy nuclear structure, and the one most directly traceable to the funda- 
mental many-body problem. Difficult though it may be, solving the shell-model 
problem is of fundamental importance to our understanding of the correlations 
found in nuclei. 

One avenue of research during the past few years has been in the area of the nu- 
clear shell model solved not by diagonalization, but by integration. In what follows, 



we will describe the shell-model Monte Carlo (SMMC) method and discuss several 
recent and interesting results obtained from theory. These include calculations in 
sd-pf '-shell neutron-rich nuclei, a discussion of electron-capture rates in /p-shell 
nuclei, pairing correlations in medium-mass nuclei near N=Z, and studies of level 
densities in rare-earth nuclei. 



II SMMC METHODS 

In the following we briefly outline the formalism of the SMMC method. We begin 
with a brief description of statistical mechanics techniques used in our approach, 
then discuss the Hubbard-Stratonovich transformation, and end with a discussion 
of Monte Carlo sampling procedures. We refer the reader to previous works [1,3] 
for a more detailed exposition. 



A Observables 

SMMC methods rely on an ability to calculate the imaginary-time many-body 
evolution operator, exp(— j3H), where (3 is a real c-number. The many-body Hamil- 
tonian can be written schematically as 

H = eO + -VOO , (1) 

where O is a density operator, V is the strength of the two-body interaction, and 
e is a single-particle energy. In the full problem, there are many such quantities 
with various orbital indices that are summed over, but we omit them here for the 
sake of clarity. 

While the SMMC technique does not result in a complete solution to the many- 
body problem in the sense of giving all eigenvalues and eigenstates of H, it can 
result in much useful information. For example, the expectation value of some 
observable, Q, can be obtained by calculating 

Tre-f> H n 

Here, j3 = T _1 is interpreted as the inverse of the temperature T, and the many- 
body trace is defined as 

TrX = £(*|X| ? >, (3) 

i 

where the sum is over many-body states of the system. In the canonical ensemble, 
this sum is over all states with a specified number of nucleons (implemented by 
number projection [2,3]), while the grand canonical ensemble introduces a chemical 
potential and sums over all many-body states. 



In the limit of low temperature (T — > or (3 — > oo), the canonical trace reduces 
to a ground-state expectation value. Alternatively, if |$) is a many-body trial state 
not orthogonal to the exact ground state, |^), then e~@ H can be used as a filter to 
refine |<&) to |^) as (3 becomes large. An observable can be calculated in this "zero 
temperature" method as 

($|e-^|$) (*|*) ' 1 ' 

If is the Hamiltonian, then (4) at (3 = is the variational estimate of the energy, 
and improves as (3 increases. Of course, the efficiency of the refinement for any 
observable depends upon the degree to which |$) approximates |\&). 

Beyond such static properties, e~^ H allows us to obtain some information about 
the dynamical response of the system. For an operator Q, the response function, 
Rq(t), in the canonical ensemble is defined as 

Rn(r) = Ti J H e = (nt(r)n(0)>, (5) 

where fit(r) = e TH ^e~ rH is the imaginary-time Heisenberg operator. Interesting 
choices for Q are the annihiliation operators for particular orbitals, the Gamow- 
Teller, Ml, or quadrupole moment, etc. Inserting complete sets of A-body eigen- 
states of H ({|i), |/)} with energies E i: f) shows that 

Rn(r)^^-" Ei \(f\m\ 2 e- r ^- Ei \ (6) 

if 

where Z = J2i e~^ El is the partition function. Thus, Rq(t) is the Laplace transform 
of the strength function Sq(E): 

/oo 
e- TE S n (E)dE ; (7) 
-oo 

Sn(E) = ^e-^\(fm\ 2 5(E-E, + E i ) . (8) 

fi 

Hence, if we can calculate ifo(r), Sq(E) can be determined. Short of a full inversion 
of the Laplace transform (which is often numerically difficult), the behavior of 
Rq(t) for small r gives information about the energy- weighted moments of Sq. In 
particular, 

/oo 1 
S n (E)dE = - ]T e"^ | (f\m\ 2 = (rfn) A (9) 
-oo Zj 



is the total strength, 



*n(0) 



S n (E)EdE = \Y.^\m\mE f - E.) 



(10) 



— oo 



if 



is the first moment (the prime denotes differentiation with respect to r). 

It is important to note that we usually cannot obtain detailed spectroscopic 
information from SMMC calculations. Rather, we can calculate expectation values 
of operators in the thermodynamic ensembles or the ground state. Occasionally, 
these can indirectly furnish properties of excited states. For example, if there is 
a collective 2 + state absorbing most of the E2 strength, then the centroid of the 
quadrupole response function will be a good estimate of its energy. But, in general, 
we are without the numerous specific excitation energies and wave functions that 
characterize a direct diagonalization. This is both a blessing and a curse. The 
former is that for the very large model spaces of interest, there is no way in which 
we can deal explicitly with all of the wave functions and excitation energies. Indeed, 
we often don't need to, as experiments only measure average nuclear properties at 
a given excitation energy. The curse is that comparison with detailed properties 
of specific levels is difficult. In this sense, the SMMC method is complementary to 
direct diagonalization for modest model spaces, but is the only method for treating 
very large problems. 

B The Hubbard- St rat onovich transformation 

It remains to describe the Hubbard-Stratonovich "trick" by which e~^ H is man- 
aged. In broad terms, the difficult many-body evolution is replaced by a superpo- 
sition of an infinity of tractable one-body evolutions, each in a different external 
field, a. Integration over the external fields then reduces the many-body problem 
to quadrature. 

To illustrate the approach, let us assume that only one operator O appears in the 
Hamiltonian (1). Then all of the difficulty arises from the two-body interaction, 
that term in H quadratic in O. If H were solely linear in O, we would have a 
one-body quantum system, which is readily dealt with. To linearize the evolution, 
we employ the Gaussian identity 



Here, h is a one-body operator associated with a c-number field a, and the many- 
body evolution is obtained by integrating the one-body evolution, U a = e _/3ft- , 
over all a with a Gaussian weight. The phase, s, is 1 if V < 0, or % if V > 0. 
Equation (11) is easily verified by completing the square in the exponent of the 
integrand; since we have assumed that there is only a single operator O, there is 
no need to worry about non-commutation. 

For a realistic Hamiltonian, there will be many non-commuting density operators, 
O a , present, but we can always reduce the two-body term to diagonal form. Thus 





for a general two-body interaction in a general time-reversal invariant form, we 
write 



H = E (<<5« + ^O a ) + \ ]T V a {O a , 6 a } , (12) 

a a 

where O a is the time reverse of O a . Since, in general, [O a , Op] ^ 0, we must split 
the interval (3 into N t "time slices" of length A(3 = f3/N t) 

e-? H = [e-^ H ] N \ (13) 



and for each time slice n = l,...,N t perform a linearization similar to Eq. 11 
using auxiliary fields a an . Note that because the various O a need not commute, 
the representation of e~ Al3h must be accurate through order (A/3) 2 to achieve an 
overall accuracy of order A{3. 

We are now able to write expressions for observables as the ratio of two field 
integrals. Thus expectations of observables can be written as 



where 



and 



{ ' JVaW a ' 1 j 



TrUM ^ [ Ap\V a \ \ 

a = ~tTu7 ; = n II d ^nda an y-^- j , (15) 

U a = U Nt ...U 2 U 1 ; f/„ = e' A ^; 

K = E ( £ a + SaK^an) O a + (e a + S a V a (J* an ) O a . (16) 



This is, of course, a discrete version of a path integral over a. Because there is a 
field variable for each operator at each time slice, the dimension of the integrals Va 
can be very large, often exceeding 10 5 . The errors in Eq. 14 are of order A/3, so that 
high accuracy requires large N t and perhaps extrapolation to Nt = oo (A(3 = 0). 

Thus, the many-body observable is the weighted average (weight W a ) of the 
observable fl a calculated in an ensemble involving only the one-body evolution U a . 
Similar expressions involving two a fields (one each for e~ rH and e^^~ T ^ H ) can 
be written down for the response function (5), and all are readily adapted to the 
canonical or grand canonical ensembles or to the zero-temperature case. 

An expression of the form (14) has a number of attractive features. First, the 
problem has been reduced to quadrature — we need only calculate the ratio of two 
integrals. Second, all of the quantum mechanics (which appears in Q a ) is of the 
one-body variety, which is simply handled by the algebra of N s x iV s matrices. The 
price to pay is that we must treat the one-body problem for all possible a fields. 



C Monte Carlo quadrature and the sign problem 



We employ the Metropolis, Rosenbluth, Rosenbluth, Teller, and Teller algorithm 
[4] to generate the field configurations, a, which requires only the ability to calculate 
the weight function for a given value of the integration variables. This method 
requires that the weight function W a must be real and non-negative. Unfortunately, 
many of the Hamiltonians of physical interest suffer from a sign problem, in that 
W a is negative over significant fractions of the integration volume. The fractional 
variance of a given expectation value becomes unacceptably large as the average 
sign approaches zero. 

It was shown [3] that for even-even and N = Z nuclei there is no sign problem 
for Hamiltonians if all V a < 0. Such forces include reasonable approximations to 
the realistic Hamiltonian like pairing-plus-multipole interactions. However, for an 
arbitrary Hamiltonian, we are not guaranteed that all V a < (see, for example, 
Alhassid et al. [6]). However, we may expect that a realistic Hamiltonian will 
be dominated by terms like those of the schematic pairing-plus-multipole force 
(which is, after all, why the schematic forces were developed) so that it is, in some 
sense, close to a Hamiltonian for which the MC is directly applicable. Thus, the 
"practical solution" to the sign problem presented in Alhassid et al. [6] is based on 
an extrapolation of observables calculated for a "nearby" family of Hamiltonians 
whose integrands have a positive sign. Success depends crucially upon the degree 
of extrapolation required. Empirically, one finds that for all of the many realistic 
interactions tested in the sd- and p /-shells, the extrapolation required is modest, 
amounting to a factor-of-two variation in the isovector monopole pairing strength. 

Based on the above observation, it is possible to decompose H in Eq. 12 into its 
"good" and "bad" parts, H = Hq + Hb- The "good" Hamiltonian, Hq, includes, 
in addition to the one-body terms, all the two-body interactions with V a < 0, while 
the "bad" Hamiltonian, H B , contains all interactions with V a > 0. By construction, 
calculations with Hq alone have = 1 and are thus free of the sign problem. 

We define a family of Hamiltonians, H g , that depend on a continuous real pa- 
rameter g as H g = f(g)Hc + gHs, so that H g =\ = H, and f(g) is a function with 
/(l) = 1 and f(g < 0) > that can be chosen to make the extrapolations less 
severe. (In practical applications, f(g) = 1 — (1 — g)/x with x ~ 4, and applied 
only to the two-body terms in Hq has been found to be a good choice.) If the V a 
that are large in magnitude are "good," we expect that H g=0 = Hq is a reasonable 
starting point for the calculation of an observable (Q). One might then hope to 
calculate (Q) g = Tr (fie _/3iffl )/Tr (e"^ 9 ) for small g > and then to extrapolate to 
g = 1, but typically (<&) collapses even for small positive g. However, it is evident 
from our construction that H g is characterized by $ CT = 1 for any g < 0, since all 
the "bad" V a (> 0) are replaced by "good" gV a < 0. We can therefore calculate 
{Q) g for any g < by a Monte Carlo sampling that is free of the sign problem. If 
{Q) g is a smooth function of g, it should then be possible to extrapolate to g — 1 
(i.e., to the original Hamiltonian) from g < 0. We emphasize that g — is not 
expected to be a singular point of (fi) 9 ; it is special only in the Monte Carlo eval- 



uation. The extrapolation methods we employ have been tested against standard 
shell-model diagonalizations in many cases, and have, in general, been shown to 
work very well [3]. 



Studies of extremely neutron-rich nuclei have revealed a number of intriguing 
new phenomena. Two sets of these nuclei that have received particular attention 
are those with neutron number N in the vicinity of the IsOd and O/7/2 shell closures 
(N ~ 20 and N « 28). Experimental studies of neutron-rich Mg and Na isotopes 
indicate the onset of deformation, as well as the modification of the N = 20 shell 
gap for 32 Mg and nearby nuclei [7]. Inspired by the rich set of phenomena occurring 
near the N = 20 shell closure when N ^> Z, attention has been directed to nuclei 
near the N = 28 (sub) shell closure for a number of S and Ar isotopes [8,9] where 
similar, but less dramatic, effects have been seen as well. 

In parallel with the experimental efforts, there have been several theoretical stud- 
ies seeking to understand and, in some cases, predict properties of these unstable 
nuclei. Both mean-field [10,11] and shell-model calculations [8,9,12-16] have been 
proposed. The latter require a severe truncation to achieve tractable model spaces, 
since the successful description of these nuclei involves active nucleons in both the 
sd- and the ^/-shells. The natural basis for the problem is therefore the full sd-pf 
space, which puts it out of reach of exact diagonalization on current hardware. 

SMMC methods offer an alternative to direct diagonalization when the bases 
become very large. Though SMMC provides limited detailed spectroscopic in- 
formation, it can predict, with good accuracy, overall nuclear properties such as 
masses, total strengths, strength distributions, and deformation — precisely those 
quantities probed by the recent experiments. It thus seems natural to apply SMMC 
methods to these unstable neutron-rich nuclei. Two questions will arise — center- 
of-mass motion and choice of the interaction — that are not exactly new, but 
demand special treatment in very large spaces. These questions were addressed in 
detail in Ref. [17]. We present a brief selection of results here. 

There is limited experimental information about the highly unstable, neutron- 
rich nuclei under consideration. In many cases only the mass, excitation energy 
of the first excited state, the B(E2) to that state, and the /3-decay rate is known, 
and not even all of this information is available in some cases. From the measured 
B(E2), an estimate of the nuclear deformation parameter, fa, has been obtained 
via the usual relation 



III APPLICATIONS 



A sd-pf nuclei 




(17) 



with R = 1.2A 1 / 3 fm and B(E2) given in e fm . 
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FIGURE 1. Top: Difference between theoretical and experimental binding energies for the 
sd-pf -shell nuclei studied in this work. Bottom: Experimental and theoretical B(E2) values. 



Much of the interest in the region stems from the unexpectedly large values 
of the deduced fa, results which suggest the onset of deformation and have led 
to speculations about the vanishing of the N = 20 and N = 28 shell gaps. The 
lowering in energy of the 2f state supports this interpretation. The most thoroughly 
studied case, and the one which most convincingly demonstrates these phenomena, 
is 32 Mg with its extremely large B(E2) = 454 ± 78e 2 fm 4 and corresponding fa = 
0.513 [7]; however, a word of caution is necessary when deciding on the basis of 
this limited information that we are in the presence of well-deformed rotors: for 
22 Mg, we would obtain fa = 0.67, even more spectacular, and for 12 C, fa = 0.8, 
well above the superdeformed bands. 

Most of the measured observables can be calculated within the SMMC frame- 
work. It is well known that in deformed nuclei the total B(E2) strength is al- 
most saturated by the 0+ — > 2± transition (typically 80% to 90% of the strength 
lies in this transition). Thus the total strength calculated by SMMC should only 
slightly overestimate the strength of the measured transition. In Fig. 1 the SMMC 
computed B(E2, total) values are compared to the experimental B(E2;0+ S — > 2+) 
values. Reasonable agreement with experimental data across the space is obtained 
when one chooses effective charges of e p = 1.5 and e n = 0.5. Using these same 
effective charges, the USD values for the B(E2, 0+ — > 2f) of the srf-shell nuclei 
32 Mg and 30 Ne are 177.1 and 143.2 e 2 fm 4 , respectively, far lower than the full sd- 
pf calculated and experimental values. All of the theoretical calculations require 
excitations to the p/-shell before reasonable values can be obtained. We note a 



TABLE 1. Comparisons of the SMMC electron 
capture rates with the total (A oc ) and partial 
Gamow- Teller (A^. T ) rates as given in Rcf. [19]. Phys- 
ical conditions at which the comparisons were made 
are p 7 = 5.86, T 9 = 3.40, and Y e = 0.47 for the up- 
per part of the table, and pi = 10.7, T 9 = 3.65, and 
Y e = 0.455 for the lower part. 



Nucleus 


A cc (sec -1 ) 


A oc (sec -1 ) 


A£ J (sec- 1 ) 




(SMMC) 


(Ref. [19]) 


(Ref. [19]) 


55 Co 


3.89E-04 


1.41E-01 


1.23E-01 


57 Co 


3.34E-06 


3.50E-03 


1.31E-04 


55p e 


1.20E-08 


1.61E-03 


1.16E-07 


56 Ni 


3.47E-02 


1.60E-02 


6.34E-03 


58 Ni 


1.01E-03 


6.36E-04 


4.04E-06 


60 Ni 


7.39E-05 


1.49E-06 


4.86E-07 


5y Co 


3.44E-07 


2.09E-04 


6.37E-05 


57 Co 


2.06E-05 


7.65E-03 


3.69E-04 


55p e 


1.07E-07 


3.80E-03 


5.51E-07 


56p e 


9.80E-06 


4.68E-07 


6.60E-10 


54p e 


3.84E-04 


9.50E-04 


3.85E-06 


51y 


1.06E-06 


1.24E-05 


9.46E-09 


52 Cr 


1.32E-04 


2.01E-07 


1.59E-10 


60 Ni 


3.61E-04 


7.64E-06 


2.12E-06 



general agreement among all calculations of the B(E2) for 46 Ar, although they are 
typically larger than experimental data would suggest. We also note a somewhat 
lower value of the B(E2) in this calculation as compared to experiment and other 
theoretical calculations in the case of 42 S [15]. 

Also shown in Fig. 1 are the differences between experimental and theoretical 
binding energies for nuclei in this region. Agreement is quite good overall. Further 
details of the interaction and results may be found in [17]. 

B Electron capture rates for Fe-region nuclei 

The impact of nuclear structure on astrophysics has become increasingly im- 
portant, particularly in the fascinating, and presently unsolved, problem of type- 
II supernovae explosions. One key ingredient of the precollapse scenario is the 
electron-capture cross section on nuclei [18,19]. An important contribution to 
electron-capture cross sections in supernovae environments is the Gamow- Teller 
(GT) strength distribution. This strength distribution, calculated in SMMC us- 
ing Eqs. (3 and 4) above, is used to find the energy-dependent cross section for 
electron capture. In order to obtain the electron-capture rates, the cross section is 
then folded with the flux of a degenerate relativistic electron gas [20]. Note that the 



Gamow- Teller distribution is calculated at the finite nuclear temperature which, in 
principle, is the same as the one for the electron gas. 

It is important to calculate the GT strength distributions reasonably accurately 
for both the total strength and the position of the main GT peak in order to have 
a quantitative estimate for the electron-capture rates. For astrophysical purposes, 
calculating the rates to within a factor of two is required. We concentrate here on 
mid-/p-shell results for the electron-capture cross sections [20]. The Kuo-Brown 
interaction [21], modified in the monopole terms by Zuker and Poves [22], was used 
throughout these p /-shell calculations. This interaction reproduces quite nicely the 
ground- and excited-state properties of mid-/p-shell nuclei [23,24], including the to- 
tal Gamow- Teller strengths and distributions, where the overall agreement between 
theory and experiment [25] is quite reasonable. The SMMC technique allows one 
to probe the complete Oku /p-shell region without any parameter adjustments to 
the Hamiltonian, although the Gamow- Teller operator has been renormalized by 
the standard factor of 0.8. 

Do the electron-capture rates presented here indicate potential implications for 
the precollapse evolution of a type II supernova? To make a judgment on this 
important question, one should compare in Table I the SMMC rates for selected 
nuclei with those currently used in collapse calculations [19]. For the comparison, 
we choose the same physical conditions as assumed in Tables 4-6 in Aufderheide et 
al. [19]. Table I also lists the partial electron-capture rate which has been attributed 
to Gamow- Teller transitions [19]. Note that for even-parent nuclei, the present rate 
approximately agrees with the currently recommended total rate. A closer inspec- 
tion, however, shows significant differences between the present rate and the one 
attributed to the Gamow- Teller transition in Aufderheide et al. [19]. The origin of 
this discrepancy is due to the fact that Fuller et al. [26] places the Gamow- Teller 
resonance for even-even nuclei systematically at too high an excitation energy. 
This shortcoming has been corrected in Fuller et al. [26] and Aufderheide et al. 
[19] by adding an experimentally known low-lying strength in addition to the one 
attributed to Gamow- Teller transitions. However, the overall good agreement be- 
tween the SMMC results for even-even nuclei and the recommended rates indicates 
that the SMMC approach also accounts correctly for this low-lying strength. This 
has already been deduced from the good agreement between SMMC Gamow- Teller 
distributions and data including the low-energy regime [25]. 

Thus, for even-even nuclei, the SMMC approach is able to predict the total 
electron-capture rate rather reliably, even if no experimental data are available. 
Note that the SMMC rate is somewhat larger than the recommended rate for 56 Fe 
and 60 Ni. In both cases, the experimental Gamow- Teller distribution is known and 
agrees well with the SMMC results [25] . While the proposed increase of the rate for 
60 Ni is not expected to have noticeable influence on the pre-collapse evolution, the 
increased rate for 56 Fe makes this nucleus an important contributor in the change 
of Y e during the collapse (see Table 15 of Aufderheide et al. [19]). 

For electron capture on odd- A nuclei, observe that the SMMC rates, derived from 
the Gamow- Teller distributions, are significantly smaller than the recommended to- 
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FIGURE 2. Left: proton pair strength distribution for 46 V. Right: proton-neutron (T=l) pair 
strength distribution. SMMC: thick line; direct diagonalization: impulses. 



tal rate. This is due to the fact that for odd- A nuclei the Gamow- Teller transition 
peaks at rather high excitation energies in the daughter nucleus. The electron- 
capture rate on odd- A nuclei is therefore carried by weak transitions at low exci- 
tation energies. Comparing the SMMC rates to those attributed to Gamow- Teller 
transitions in Fuller et al. [26] and Aufderheide et al. [19] reveals that the latter 
have been, in general, significantly overestimated, which is caused by the fact that 
the position of the Gamow- Teller resonance is usually put at too low excitation 
energies in the daughter. The SMMC calculation implies that the Gamow- Teller 
transitions should not contribute noticeably to the electron-capture rates on odd- A 
nuclei at the low temperatures studied in Tables 14-16 in Aufderheide et al. [19]. 
Thus, the rates for odd- A nuclei given in these tables should generally be replaced 
by the non-Gamow- Teller fraction. 



C Pair correlations in nuclei 

We would now like to turn to the subject of pair correlations in nuclei, and 
calculations aimed at their understanding. Nuclei near N=Z offer a unique place 
to study proton- neutron pairing, particularly in the isospin T=l channel. In fact, 
most heavy odd-odd, N=Z nuclei beyond 40 Ca have total spin J=0, T=l ground 
states. Theoretical studies have shown that many of these nuclei have enhanced 
T=l proton- neutron correlations when compared to their even-even counterparts. 
These correlations are, to a lesser extent, present in even-even systems, but tend 
to decrease as one moves away from N=Z. In at least one nucleus in the mass 70 
region, 74 Rb, there is experimental evidence for a ground state T = 1 band [27]. 

Experimentally, pair correlations can best be measured by pair transfer on nu- 
clei. Although total cross sections are typically underpredicted when one employs 
spectroscopic factors computed from the shell model, relative two-nucleon spectro- 



scopic factors within one nucleus are more reliable. Therefore, it is necessary for 
one to calculate and measure pair transfer from both the ground and excited states 
in a nucleus. 

The SMMC method may be used to calculate the strength distribution of the pair 
annihilation operator Ajtt z , as defined in Koonin et al. [3]. The total strength of 
these pairing operators, i.e. the expectation {Aj T A jt), has been studied previously 
as a function of mass, temperature [28,29], and rotation [30]. We would like to 
briefly present here the strength distributions of the pair operators as calculated 
in SMMC. The strength distribution for the pair transfer spectroscopic factors is 
proportional to (A — 2 | A JT \ A) and is calculated by the inversion of Eq. (4). 

In future work, we will discuss the strength distributions in detail. Here we would 
like to briefly conclude by demonstrating that the SMMC results and the direct 
diagonalization results agree very nicely for the proton pair strength distributions 
in the ground state of 46 V. This is demonstrated in the left panel of Fig. 2. Shown 
in the right panel is the isovector proton-neutron pairing strength distribution with 
respect to the daughter nucleus. Notice that the overall strength is much larger 
in the proton-neutron channel, as discussed previously in Langanke et al. [29], and 
that the peak is several MeV lower in excitation relative to the like-particle channel. 
In both cases the strength distribution in 46 V differs significantly from that found 
in 48 Cr, where one finds that the dominant component is a ground-state to ground- 
state transition involving mainly particles in the O/7/2 single-particle state. In both 
odd-odd N=Z channels, the distribution is fairly highly fragmented. 

D Rare earth nuclei 

We have recently applied SMMC techniques to survey rare-earth nuclei in the Dy 
region. This extensive study formed the thesis topic of J. A. White [36], whose goal 
was to examine how the phenomenologically motivated "pairing-plus-quadrupole" 
interaction compares in exact shell-model solutions with other methods. We also 
examined how the shell-model solutions compare with experimental data; static 
path approximation (SPA) calculations were also performed. There have been ef- 
forts recently by others to use SPA calculations, since they are simpler and faster 
(see [37,38] as examples). However, we found that SPA results are not consistently 
good. This study was also designed to investigate whether the phenomenological 
pairing-plus-quadrupole-type interactions can be used in exact solutions for large 
model spaces, and whether the interaction parameters require significant renormal- 
ization when using SPA. 

We discuss here one particular aspect of that work, namely level density calcu- 
lations. Details may be found in [39]. We used the Kumar-Baranger Hamiltonian 
with parameters appropriate for this region. Our single-particle space included the 
50-82 subshell for the protons and the 82-126 shell for the neutrons. While several 
interesting aspects of these systems were studied in SMMC, we limit our discussion 
here to the level densities obtained for 162 Dy. 



SMMC is an excellent way to calculate level densities. E((3) is calculated for 
many values of f3 which determine the partition function, Z, as 



ln[Z(/?)/Z(0)] = - [ d/3'E(/3') 
Jo 



(18) 



Z(0) is the total number of available states in the space. The level density is then 
computed as an inverse Laplace transform of Z. Here, the last step is performed 
with a saddle point approximation with j3~ 2 C = —dE/d(3: 



SMMC has been used recently to calculate level density in iron region nuclei [40], 
and here we demonstrate its use in the rare-earth region. 

The comparison of SMMC density in 162 Dy with the Tveter et al. [41] data is 
displayed in Fig. 3. The experimental method can reveal fine structure, but does 
not determine the absolute density magnitude. The SMMC calculation is scaled by 
a factor to facilitate comparison. In this case, the factor has been chosen to make 
the curves agree at lower excitation energies. From 1-3 MeV, the agreement is very 
good. From 3-5 MeV, the SMMC density increases more rapidly than the data. 
This deviation from the data cannot be accounted for by statistical errors in either 
the calculation or measurement. Near 6 MeV, the measured density briefly flattens 
before increasing and this also appears in the calculation, but the measurement 
errors are larger at that point. 

The measured density includes all states included in the theoretical calculation 
plus some others, so that one would expect the measured density to be greater 
than or equal to the calculated density and never smaller. We may have instead 
chosen our constant to match the densities for moderate excitations and let the 
measured density be higher than the SMMC density for lower energies (1-3 MeV). 
Comparing structure between SMMC and data is difficult for the lowest energies 
due to statistical errors in the calculation and comparison at the upper range of the 
SMMC calculation, i.e., E 15 MeV, is unfortunately impossible since the data 
only extend to about 8 MeV excitation energy. 



In these Proceedings, we have used four specific examples (there are several 
others) for which the SMMC calculations have proven very useful in understanding 
the properties of nuclei in systems where the number of valence particles prohibits 
the use of more traditional approaches. The method has proven to be a valuable tool 
for furthering our understanding of nuclear structure and astrophysics. Continued 
developments in both creating useful interactions and shell-model technology should 
continue to enhance our ability to understand nuclei far from stability in the coming 
years. 



S(E) =(3E + \nZ((3) 

p{E) = (27r/r 2 C)- 1/2 exp(S) 



(19) 
(20) 



IV CONCLUSIONS 
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FIGURE 3. SMMC density vs. experimental data in 162 Dy. 
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